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Abstract 

We calculate three-point functions of two protected operators and one twist-two operator 
with arbitrary even spin j in J\f — 4 super Yang-Mills theory to one-loop order. In order 
to carry out the calculations we project the indices of the spin j operator to the light-cone 
and evaluate the correlator in a soft-limit where the momentum coming in at the spin j 
operator becomes zero. This limit largely simplifies the perturbative calculation, since all 
three-point diagrams effectively reduce to two-point diagrams and the dependence on the 
one-loop mixing matrix drops out completely. The results of our direct calculation are in 
agreement with the structure constants obtained by F.A. Dolan and H. Osborn from the 
operator product expansion of four-point functions of half-BPS operators. 



1 Introduction and Conclusions 



The full solution of M = 4 SYM theory is equivalent to the ability of calculating any 
observable of interest in this quantum field theory exactly or perturbatively to arbitrary 
order in the coupling constant. In recent years, such impressive achievements have indeed 
been made for certain observables in planar M = 4 SYM, largely due to the existence of 
integrability [U El El E] in the AdS/CFT correspondence 0, El C], see [S] for a review. 

Using the methods of integrability, all-order conjectures, e.g. for the cusp anomalous 
dimension of twist two operators could be obtained [HI EI]- Using this, also for scattering 
amplitudes of four and five external particles in planar M = 4 SYM, all-order results 
became available [HI [T2] . These all-order results were achieved using the duality between 
MHV amplitudes and light-like Wilson loops [T31 El EE] , see [El HZ] for reviews. 

Given these results, it is natural to ask, whether similar advances can be made for 
three-point functions of gauge invariant operators. Not only are three-point functions 
the next natural correlators to address after two-point functions, but they possess a par- 
ticular importance due to the operator product expansion (OPE). Given the knowledge 
of all structure constants that appear in the theory, one can in principle construct any 
higher-point correlation function using the OPE. Therefore, all-order results for struc- 
ture constants, together with all-order results for anomalous dimensions yield important 
information on any higher-point function. 

The form of three-point functions is fixed by conformal symmetry and the only depen- 
dence on the coupling constant is contained in the anomalous dimensions of the operators 
and the structure constants which receive radiative corrections in the coupling 

Ca^ig 2 ) = C a $ 7 + g 2 C^ + ... (1) 

Similarly as for two-point functions, there are non-renormalisation theorems for three- 
point correlators of half-BPS operators, which guarantee that they do not get quantum 
corrections [TBI US 1201 ED E21 E3] • 

Direct computations of three-point functions in M = 4 SYM theory have been per- 
formed injSaESlESlEZlEHlESlEniEIlES]. In [33] a large number of three-point functions 
involving scalar primary operators of up to and including length five is considered. 

The role of integrability for three-point function calculations was first addressed in [2S1 
EZ1ES]- Applications of integrability methods can be found in [511 E3 ESI EZ] at tree- level 
and in [3E1 EZ] for loop-level three-point functions of scalar single trace operators. It does 
indeed turn out, that three-point functions can be studied efficiently using integrability. 

The question of a weak-strong coupling matching was addressed in [23 ESI HD1 El H21 
145] . The SL(2) sector was adressed in [43] and agreement between the structure constants 
at weak and strong coupling was found, where a BPS state and two operators with large 
spin were considered at tree-level. It would be interesting to promote this study to the 
loop-level, which requires the knowledge of the correlator of one BPS operator and two 
operators with spin at one-loop level. Furthermore, it would be interesting to see, how one 
can make use of integrability methods in order to calculate this correlator at loop-level. 

Here, we make a first step towards this direction by providing a direct field theory 
computation of three-point correlators with two protected operators and one twist-two 
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operator. This calculation serves as a starting point for the case of two operators with 
spin, where several distinct space-time structures enter with a priori different structure 
constants, see e.g. [3S]- A direct calculation is possibly more suitable for the devel- 
opment of integrability methods than the extraction of the same coefficient from the OPE. 

More explicitly, here we calculate three-point functions in M = 4 SYM theory involv- 
ing two protected scalar operators 

0(x) = Tr (fa(x)fa 3 (x)) , 6(x) = Tr ^ 12 (x)<p 13 (x)) (2) 

and one twist-two operator with arbitrary even spin j, which is totally symmetric and 
traceless in all indices and schematically of the form 

0^.(x) = Tr (d^.D^J 12 (x)D^ +1 ..D^ 12 (x)) + ... (3) 

where the ellipses stand for a certain distribution of the derivatives, which is given explic- 
itly in section EJ Conformal symmetry fixes the three-point functions up to the structure 
constants C ^. We use the light-cone projection, see appendix IA.21 where all indices are 
contracted with a light-like vector z 11 . Then the correlator with the renormalised operator 
see (EI), reads 
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where 6 = A,- — j is the twist (dimension minus spin) of the operator O 



13 _ f_23 = 1 o M i / f_2 
2 2 n u x 3 I 2 

13 x 23 z \ X 13, 



*W = r% , ^(x 13 , x 23 ) = ^-^ = > X3 ln[^\ (5) 



and where xfj = xf — Xj. Note, that the dimension of the twist operator and the structure 
constants both acquire corrections in perturbation theory 

9 = ^-3 = 2 + 7j (s 2 ), C o6j (g*) = C®. + g 2 C% 3 + Otf) . (6) 

In order to carry out the calculations we evaluate the correlator in a soft-limit, where 
the momentum coming in at the spin j operator becomes zero. In position space this 
corresponds to an integration over the corresponding point X3 

''" P e^(0{p)d{-p)Q 3 m = ( d d x 3 (0(x 1 )6(x 2 )6 j (x 3 )) . (7) 



(27T) 



This limit largely simplifies the perturbative calculation, since all three-point diagrams 
effectively reduce to two-point diagrams and the dependence on the one-loop mixing 
matrix drops out completely, as we will see in section EJ We find that the correction to 
the structure constant takes the simple form 



C'oojig) = C% (l + §^ miHj - H 2j ) - H h2 ) + OU, 
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where Hj >m are generalised harmonic sums Hj = Y? n V n > Hj,r — J2 J n ^/n r . The result is in 
agreement with the extraction of the structure constants from the analysis of the operator 
product expansion (OPE) of four-point functions of half-BPS operators [IB]. Recently, 
the OPE and the three-loop expression for the four-point correlator of half-BPS operators 
|4"T] were used to extract the structure constants up to three loops [IB]. Our calculation 
thus directly confirms the results obtained from the OPE. 

In the following sections we will only summarise the main steps of our calculation in 
order to keep the presentation short, more details will be made available in [49J. 



2 Three-Point Functions of Twist Operators 

We use the mostly minus metric and the conventions as given in Appendix |A] and adapt 
the notation from [30] ■ The renormalised twist-two operators are given by 

6 j (x) = J Ez jk &- k 6 k (x), (9) 

k 

where the tree-level correlator is 

= E *)k Tr (D k <p 12 (x)&- k <p 12 (*)) , (10) 

fc=0 

1 /2 _____ 

the coefficients a- k are related to the Gegenbauer polynomials as in (1571) and the one-loop 
mixing matrix in the conformal scheme has the form 

Z jk = 8 jk + g 2 zf k ] + O(g') = 8 jk + g 2 + U jk zf^ + O(g') . (11) 

At one-loop level the divergent part of the renormalisation matrix is diagonal with Zj = 
Hj/{A r K d l 2 ) exp(e7s) and determines the anomalous dimension via jj = —fid/dfilnZj(/j,). 
The finite matrix Bj k (g 2 ) diagonalizes the anomalous dimension matrix at higher- loop 
level '~Yj{g 2 )5jk = (B~ 1/ -fB)j k and accounts for non-diagonal contributions to the two-point 
functions at loop-level. We do not specify Bj k , since it drops out in the limit that we 
consider. The operators given by ([§]) then have diagonal, conformal two-point functions 
also at loop- level and are fixed up to the normalization Cj(g 2 ) 

(Q^OQ^)) = 5 jk Cj(g 2 ) (j even) , (12) 

where 9 = Aj — j = 2 + , jj(g 2 ) and 1^ = ~ 2x^2^12/^12- m order to read-off the 
structure constants in the limit (j7J we have to integrate (jlj) over x 3 , which yields 

/ d^OOO,) = N{g 2 ) (Cg> . + g 2 C^) ^^,2 ( 13 ) 
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where N(g 2 ) is a normalisation factor explicitly given in (I45p . By calculating the left-hand 
side in the limit pi + P2 — > one can thus easily read-off the structure constants, after 
Fourier transforming the momentum space expression^ 



3 Tree-Level Calculation 

In momentum space, the tree-level three-point function with the twist-operator in the 
representation (1551) reads 

1 r d d k (Pi +P2) j C 1/2 ( 2k 7 pl 7 h ) 

Now we take the limit p% + p 2 — > in momentum space. Then, due to the factor (pi 
only the term with the highest power, i.e. j, in the Gegenbauer polynomial can survive. 
The corresponding coefficient reads 



# = ^FTTWXTT cy\ x) = ±c%V. (15) 

fc=0 



I )!'(./ 



Thus, in this limit the three-point integral in f)14p becomes a two-point integral with a 
doubled propagator 

d*k (Pl+P2)< 2 (^^) - 1/2/ - <ffc (fcV 

(27r) d A; 2 (j9 1 -fc) 2 (p 1 +p 2 -A;) 2 JJ 7 (2ir) d k\ Pl - kf ' 1 J 

which can easily be solved using (1421) in the appendix. Fourier transformation to x-space 
using (14 1|) and insertion of the coefficients yields 

d d x (odd ) - ia e 5aa ^- 7+j mnh-in 3 -2 + h) 

Comparing with (1131) we can directly read off the tree-level structure constant for d = 4 
and arbitrary even spin j 

C oo 3 = 9 6Sm2 ^ r (j + 1) ■ (18) 



4 One-Loop Calculation 

In order to read-off the one-loop structure constant we have to compute the renormalised 
three-point function using the renormalised operator given by ([§l)- (fTTj) . From 

/ d*x z (066 j ) = £% fc / d*x z %- k (OOd k ) =Y,Z jk 5 jk [ d d x 3 (OOd k ) 

J k J k 

1 An analogue of this procedure appears in [5T] , where a deformation of the N = 4 SYM Lagrangian 
by integrated local operators is considered and two-point functions of the deformed theory are studied at 
leading order in the deformation parameter. 
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it is however immediately clear, that the finite matrix Bjf. drops out in the limit that 
we consider, which vastly simplifies the calculation. As can be seen from fllip . we then 
only need to calculate the one-loop diagrams and subtract the tree-level contribution^! 
multiplied with the renormalisation constant Zj. The one-loop diagrams that enter the 
calculation are shown in figure [TJ 
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Figure 1: Feynman diagrams contributing to the three-point function at one loop. For j 
even the diagrams in the second and third row are identical. 



It is instructive to write down all diagrams before employing the limit p\ +P2 — > and 
to convince oneself that diagrams which could lead to ambiguities, such as the diagram 
in figure [Tgl which vanishes in dimensional regularisation when taking the limit, cancel 
with contributions from other diagrams. 

Indeed, all diagrams with a scalar four-point interaction vertex, i.e. fig. [Tel [Tg| [Ik 
cancel against contributions from fig. [Tbj [Tfj [Tjj which becomes clear when considering 



the decomposition of diagrams with two gluon vertices into scalar integrals, see figure El 




-2pl 








Figure 2: Due to the momenta from the gluon vertices, the integrand decomposes into 
simpler integrals, which cancel with the self-energy and four-scalar interaction terms. 



For j = 0, the diagrams in the last column do not appear, since there is no covariant 
derivative that the gauge field at X3 could arise from and half of the diagrams [Tel E3 
cancels against a contribution from [TfJ, half of the diagrams [Taj [TT] cancels against a 
contribution from [TTJ and half of the diagrams [Tel HJJ cancels against a contribution from 
[Tbl The remaining finite contributions from [Tbl [If] [Tj| cancel due to a relation between 



2 It needs to be evaluated including 0(e) terms, because it is multiplied with 1/e. 
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scalar three-point integrals given in [52]. Thus for j = all contributions exactly cancel, 
which is due to the fact that in this case all of the operators are protected. 

For j ^ the above cancellations between divergent diagrams where the divergence 
is located at the points x\ resp. x 2 , i.e. where the BPS operators sit, remain true as one 
would expectH. The cancellations between half of the diagrams [Tel M with a contribution 
from [lb] does however not take place anymore, and the spin-j operator at 23 acquires an 
anomalous dimension jj. For the remaining diagrams in the limit p\ + P2, for the same 

1 /2 

reason as in the tree-level calculation, we only need the coefficient with the highest 
power of the Gegenbauer polynomial that was given in 0151) . Since pi = —p2 = p and j is 
even we find the same result for all diagrams which are equal under p\ -H- P2 ■ The integrals 
turn into simple two-point integrals that can be easily solved using (I42p . Denoting from 
now on pi = p we find 

(oo6 j f^+(oo6 j f3= (19) 



-c 



'33 y "3\^ 2' ;V JV ' y 3\ ' " (—p2\h-d " 

where b n (ai, 02) is defined in ( 143ft . As mentioned before, half of the self-energy diagrams 
in fig. [Tel El do not cancel for j^O and the expression is 

hood j j^+ho66 j 'P = (20) 



4 2 i^g*N5™V-%(l, 1)6,(4 - ^ l)^_ d 



Applying the limit to the diagrams where the divergence is located at X\ we find 

(006^+ (006^+ hpod 3 f^+ -(oddjf^ (21) 

= i^gtNS^- 2 ^ 2 ^ 



33 (_p2\d~5 



(2 % (1, 1, 1, 2, 1) + % (1, 1, 1, 1, 1) + % (1, 1, 1, 2, 0)) . 

The integrals c nm (ai, a 2 , 03, 04, as) are defined in (T4T1) and solvecQ in Appendix ICl using 
the IBP technique. All these integrals are finite. The diagrams where the divergence is 
located at x 2 , i.e. [Ik], [TJJ da] [TT] yield the same result, since they are related to the ones 
above by pi <->j) 2 . 

For j ^ we have additional diagrams with exactly one gauge field from the covariant 
derivative in Oj. The diagrams are shown in figure [Td| [Th| [TTl Diagrams [lh] and [TTJ are 



3 Would divergences at these points remain, the operators acquired an anomalous dimension, we do 
however have protected operators with vanishing anomalous dimension at x\ , X2 ■ 
4 These integrals were considered before in [53] , 
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equal in our limit. The appearing integrals are simple two-point integrals and we find 

(oMdhMoriBL -i«-^M4 - d/2, i) ± # ± £ (*) ( m ; l ) (-1: 

4 fc=l m =l n=0 \ m / \ 71 / 

( _^ )5 - d 1) - h-n-l(h 1) + 6i-fc+n(l, 1) ~ &j-fc+n+l(l, 1)] • (22) 

and one can proceed to solve the sums. Diagram [id] is calculated in a very similar way 
and reads 

Z k=l m=l n=0 V 71 / \ 71 / 

p* 



:-p 2 ) 5 - 



^Ck—l—nj—k+n C-k—n,j—k+n Qc— 1— n,j— k+n+l ) j (23) 



where all integrals c nm = c nm (l, 1, 1, 1, 1) are finite and the solution can be found in 
appendix O 

4.1 Full Bare Three-Point Function 

Taking into account the exact cancellations between diagrams, as well as the fact that 
diagrams with p\ -H- p2 are identical in the limit p\ — > —p2 for j even, the bare one-loop 
contribution to the three-point function is given by 

(OOdj)^ = (OOOjp 1 (24) 

a=a..l 

= (H25 + (120]) + 2 x {JZID + 2 x + dZ3D • 

We have calculated all diagrams in momentum space. In order to read off the structure 
constant we Fourier transform the expression to position space using (|4ip . 

As a check of the calculation one can extract the divergent part of (12~4"1) and read 
off the anomalous dimension from the three-point function. As mentioned before, only 
divergences located at x 3 remain and we get the following contributions to the anomalous 
dimension 

JIIWB/2 = ^V mm_ g 2 N ( 1_\ [T7IMU1 9 2 N ( 1 \ 

lj ^ 4tt 2 ' ' 4tt 2 ^ j + \) ' ~ 4tt 2 \ 3 j + 1 J ' 

such that we recover the well-known one-loop contribution to the anomalous dimension 
of twist-two operators 



7j = [^J *Hs + 0{g A ) • (25) 

Using the tree-level three-point function and (jHj), ffTTl) . and (1451) we find that the 
one-loop correction to the structure constant is 

C ooJ C o6 3 = ( 5 ^) 2 " WWW) - f ^) , (26) 
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(a) 





(c) (d) (e) 

Figure 3: One-loop corrections to the two-point function of (OjOj). 



4.2 Normalisation Invariant Structure Constants 

For the structure constants appearing in the operator product expansion it is useful to 
take the normalisation of the two-point functions to be one. The operators (I38p in terms 
of Gegenbauer polynomials are not normalised to one, but they have the perturbative 
expansion 



i (x 1 )d,(x 2 ))=^ (Cf +g 2 C? ) )2V 



+ 0(g 4 



(27) 



and one can explicitly calculate C\ , C- . Using the Schwinger parametrisation of the 



propagator and some properties of the Gegenbauer polynomials as e.g. in 
that the operators (I38p have the tree-level normalisation 



(0) 



4 r(2j + i) 



4 r 

9 o 



2 5 7T 4 



one finds 



(2f 



At one-loop level the diagrams shown in figure [3] appear. The calculation is well- 
known, technically similar to the three-point calculation and all appearing integrals can 
be found in the appendix We shall not reproduce the complete calculation here, but 
refer the reader to [3H] for more details on this and all other calculations. For the one-loop 
result of the normalisation we find 

9 2 N 



Jj " 4vr 2 



(3H(j f - 2H(j)H(2j)) 



(29) 



The BPS operators ([2]) have the normalisation C = 2 6 /7T 4 and do not get quantum 
corrections. We can define a normalisation invariant structure constant CL,?,. which is 

OOj 

related to C q, through 



^OOj 



(0) 

QQj 



1+9 



2 | OOj 



1C 



OOj ~ y oOj ' 



,2^/(1) 



(30) 



and corresponds to the structure constants, where all operators are normalised to one 
Inserting f[2"9l and fl2"E|) we thus find 

9 2 N 



C'ooi(9 2 ) 



'OOjW I ~ OOj 

which is the result quoted in the introduction. 



(31) 
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A Notation and Conventions 

Throughout this paper we use the Minkowski-space metric with signature (H ...). 

A.l M = 4 SYM Lagrangian 

We use the following action of Af = 4 SYM theory 

S = ^Jd d x(- \F; v F a ^ + -D^ AB D^\ B (32) 
- Y Q f abc f ade <p b - AB ^ CD ^ e cD + f\l A {a^D^\f 

* rabc\a,aA7b \c,B , * rabc\a A.b,AB\C,a\ 
~ —^J A <PAB A a + —^J A A,a ( P A B J > 

where we have taken the trace of the matrix valued fields, e.g. $ab — 'PabT' 1 using 
the normalisation of the generators T a in the fundamental representation according to 
Tr (T a T b ) = 5 ab /2. 

A. 2 Light-Cone Projection 

The calculations can be considerably simplified by projecting the indices to the light-cone 
using a light-like vector 

6, CV..,,..:" z 2 = 0. (33) 

The indices can be recovered by repeated application of a second order differential operator 
in z M 

A„ = ((h - 1) + z ■ d)d z ^ - l -z^d z ■ d z (34) 
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in the presence of the constraint z 2 = 0, see e.g. [21]. In the following, we denote with a 
hat the contraction with z^, e.g. 

d 

x = x^z^ , d x = z^d x ^ = z^" q ^ • (35) 
The twist-two operator Oj projected to the light-cone then reads 

Oj(x) = E «J? Tr (^ fc 12 (x)^- fc 12 (^)) , (36) 

fe=0 

where D = D^z^ is the light-cone projected covariant derivative and the numerical coef- 
ficients related to the so-called G eg enbauer polynomials Cj(x) such that 

t;^V- fc = (x + i/yq(^) (37) 

and v — h — 3/2, where h — d/2. Therefore, we can rewrite the operators in the bi-local 
form [SSI EE] 



or = (9 a + B b y cf ( §^4 ) Tr (0 l2 (^)0 l2 (^) 

\O a + Ob I 



(38) 



— 



Explicitly, the coefficients are given by aj{ 2 = (— l) fc (jQ(^. These operators have confor- 
mal two-point functions which are of the form 



(OjMOkW) = 5 jk C^ f ^ 12 ' (j even) , (39) 



(-40 

and the three-point functions read 



(*i2, 3 ; 



(owoMOjM) = c odJ i2|Al+Aa _ fl . l ;;;::L 2|a;23l A 2+e -A 1 . (40) 



where \x^\ = (-a;|) 1/2 and Y 12 - 3 = Y^(x 13 ,x 23 )z fl . 

B Details of the Calculation 

B.l Fourier Transformation and Bubble Integrals 

In Minkowski-space with signature (H ...) we have 

d d p e~ ip - x r(| — k) 1 z 



(41) 



11 



The two-point integral with momenta in the numerator is 



/d^ k k 
(2iry(-k*-ie)«i(-(p + k)*-ie)<* (42) 

fin 

= b n (ai,a 2 ) 



- p 2 _ ^oi+oa-l ' 



where 



t .i)n £(j + n - ai )T(l - q 2 ) £(ai + a 2 - f ) 

6 n ( ai ,a 2 )=i—- — — f—- * — 2 . 43 

(An) 11 ' 2 T(d + n — ati — a 2 ) r(ai)r(a 2 ) 

Note, that analogous expression for the bubble integral and the Fourier transformation in 
x-space differ by minus signs, due to the different sign of ie in the propagators. 

B.2 Normalization Factor of the x<$ Integration 

Integration of fll]) over X3 yields 

y 

J OOj) ( _ x 2 2 y+d-3+ 7j (9 2 )/2 



d d x,(00® 3 ) = iY(, 2 ) (Cf 6 . + , 2 ^%) rg^p (44) 



where 



iW, d) = _ j r(g-d/2 + i)r((d- g )/2) a r( 7 + (g -i)/2)^ 

r(d-e)rG + f r(j + 0-i) 



C Integrals by using the IBP Method 

We define the following set of integrals^ 



f d d k,l (k) m {l) n 

tmn{ai, 02, 03, O4, 5 ) - J jfe 2a 1 ( p + fc )2a 2 ( Z _ k )^l^( P + l)^ ' ( } 

where k = z^k^ is the contraction with a light-like vector z 2 = 0. The integral is symmetric 
under the simultaneous exchange (771,01,02) -H- (72,04,05). In particular in this work we 
need f mn (l, 1, 1, 1, 1) and fjo(2, 1, 1, 1, 1). Since the momentum dependence is the same 
everywhere we define c mn by stripping off the equal coefficients 

/mn( Q i) a 2-, 03, a 4, O5) = c TOn (oi, a 2 , O3, a 4 , 05)- — ^ „._w • 

(47) 



3 These integrals were considered before in [53] , 
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C.l Bubble Integrals 

The integrals (1461) with one argument set to zero are easy to solve and are the building 
blocks of f mn (l, 1, 1, 1, 1) after using the IBP method, see e.g. [57] for an introduction to 
the IBP technique. We also give more details on this technique in jlH]. We have 



c mn (ai,a 2 ,0, a 4 ,a 5 ) = (-l)^< ai & m (ai, a 2 )b n (a 4 , a 5 ) (48) 

and 

d 

c mn (ai,a 2 ,a 3 , 0,4, 0) = (-l) n (-l)^ ai 6 n (a 4 , a 3 )& m+n (ai + a 3 + a 4 - 2>° 2 ) ■ 

By the symmetry of the integral fT^BI) (m, ai, a 2 ) -h- (n, a 4 , as) this also yields 

c m „(ai, 0, a 3 , a 4 , a 5 ) = c nm (a 4 , a 5 , a 3 , ai, 0) . (49) 

Furthermore, we have 

c mn (ai,a 2 ,a 3 ,0,a 5 ) = (-1)5> £E (?) f " T J )(-1) B 

j=0 i=o V-V \ 1 J 

b n ^j(a 5 , as)b m+n -i-.j(ai, a 2 + a 3 + a 5 - ^) (50) 

and c mn (0, a 2 , a 3 , a 4 , as) is related to this integral by the symmetry (m, a\, a 2 ) <H> (n, a 4 , as). 
We can perform the sum over i and find 

I n ^ _ - m)r(ai + a 2 + a 3 + a 5 - rf) 

Cmn(Ql) a 2, ^3) U, a 5 J — — — — ■ ■ ■ r (DiJ 

1 (ajjl [a.\ + a 2 + a 3 + as — a — mj 

n AA _ d 

E ■ (-l) m+n J bj{a 5 , a 3 )bj{a 2 + a 3 + a 5 - -, a x - m) . 

The sum can be performed, e.g. using Mathematica, in terms of generalised hypergeo- 
metric functions. 

C.2 Integrals f mn (l, 1, 1, 1, 1) 

First we would like to solve the integral f mn (l, 1, 1, 1, 1), given through the definition 

f ( f fkj_ (k) m (ir 

7^01,02,03,04,05; -j (27r)M k2ai{p + k)2a2{k _ l)2aH2a4{p + l)2a5 ■ w 



imn f ) 

The integral is symmetric under simultaneous exchange of {m, ai,a 2 } -H- {n,a 4 ,as} and 
can be solved by using the IBP method. 

(d + m- 4)/ m „ = bubbles(m, n) + m f m -i, n +i , (53) 
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where we abbreviate the known integrals as 

bubbles(m, n) = f mn (2, 1, 0, 1, 1) - / mn (2, 1, 1, 0, 1) (54) 
+ / mn (l,2,0,l,l)-/ mn (l,2,l,l,0). 

We can actually solve for all integrals by writing down (|53[) for m = n + 1 and n — m — 1, 
i.e. 

(d + (n + 1) - 4)/ n+ i im _i = bubbles(n + 1, m - 1) + (n + 1) / nm (55) 

and using / mn (l, 1, 1, 1, 1) = / nm (l, 1, 1, 1, 1)- Solving (|55D for /„ + i, m _i and inserting it 
into 0531) we get f nm in terms of known integrals: 

, . (n + d — 3)bubbles(m, n) + mbubblesfn + 1, m — 1) , . 

/ mn 1,1,1,1,1 = — — — — — — ■ — . 56 

(a — 4) [a — 3 + m + n) 

Note, that the coefficient on the left-hand side is of order e in d = 4 — 2e. All appearing 
integrals were solved in section IC.ll 

C.3 Integrals / mn (2, 1, 1, 1, 1) 

For the three-point function calculation we need the integral 

f h (2, i, i, i, i) - j (27 r)2rf k\p + kf{i-k)n\ P + iy ' (57) 

Using the IBP method, just as in section ICT2l for / mn (l, 1, 1, 1, 1), we find 

(d - 5 + m)/ m>ri (2, 1, 1, 1, 1) = m/ m _i )n+1 (2, 1, 1, 1, 1) + bubbles 2 (m, n) , (58) 

where 

bubbles 2 (m, n) = 2f mn (3, 1, 0, 1, 1) - 2/ mn (3, 1, 1, 0, 1) (59) 
+ / m „(2,2,0,l,l)-/ m „(2,2,l,l,0) 

and all these integrals are solved in the following section. Thus we can obtain f mn recur- 
sively from 
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